\(\int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx\) [62]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [A] (verified)
   Fricas [A] (verification not implemented)
   Sympy [B] (verification not implemented)
   Maxima [B] (verification not implemented)
   Giac [A] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 36, antiderivative size = 162 \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=\frac {5 (2 A-5 B) c^3 x}{2 a^2}+\frac {5 (2 A-5 B) c^3 \cos (e+f x)}{2 a^2 f}-\frac {a^3 (A-B) c^3 \cos ^7(e+f x)}{3 f (a+a \sin (e+f x))^5}+\frac {2 a (2 A-5 B) c^3 \cos ^5(e+f x)}{3 f (a+a \sin (e+f x))^3}+\frac {5 (2 A-5 B) c^3 \cos ^3(e+f x)}{6 f \left (a^2+a^2 \sin (e+f x)\right )} \]

[Out]

5/2*(2*A-5*B)*c^3*x/a^2+5/2*(2*A-5*B)*c^3*cos(f*x+e)/a^2/f-1/3*a^3*(A-B)*c^3*cos(f*x+e)^7/f/(a+a*sin(f*x+e))^5
+2/3*a*(2*A-5*B)*c^3*cos(f*x+e)^5/f/(a+a*sin(f*x+e))^3+5/6*(2*A-5*B)*c^3*cos(f*x+e)^3/f/(a^2+a^2*sin(f*x+e))

Rubi [A] (verified)

Time = 0.23 (sec) , antiderivative size = 162, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.167, Rules used = {3046, 2938, 2759, 2758, 2761, 8} \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=-\frac {a^3 c^3 (A-B) \cos ^7(e+f x)}{3 f (a \sin (e+f x)+a)^5}+\frac {5 c^3 (2 A-5 B) \cos (e+f x)}{2 a^2 f}+\frac {5 c^3 (2 A-5 B) \cos ^3(e+f x)}{6 f \left (a^2 \sin (e+f x)+a^2\right )}+\frac {5 c^3 x (2 A-5 B)}{2 a^2}+\frac {2 a c^3 (2 A-5 B) \cos ^5(e+f x)}{3 f (a \sin (e+f x)+a)^3} \]

[In]

Int[((A + B*Sin[e + f*x])*(c - c*Sin[e + f*x])^3)/(a + a*Sin[e + f*x])^2,x]

[Out]

(5*(2*A - 5*B)*c^3*x)/(2*a^2) + (5*(2*A - 5*B)*c^3*Cos[e + f*x])/(2*a^2*f) - (a^3*(A - B)*c^3*Cos[e + f*x]^7)/
(3*f*(a + a*Sin[e + f*x])^5) + (2*a*(2*A - 5*B)*c^3*Cos[e + f*x]^5)/(3*f*(a + a*Sin[e + f*x])^3) + (5*(2*A - 5
*B)*c^3*Cos[e + f*x]^3)/(6*f*(a^2 + a^2*Sin[e + f*x]))

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 2758

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Simp[g*(g*C
os[e + f*x])^(p - 1)*((a + b*Sin[e + f*x])^(m + 1)/(b*f*(m + p))), x] + Dist[g^2*((p - 1)/(a*(m + p))), Int[(g
*Cos[e + f*x])^(p - 2)*(a + b*Sin[e + f*x])^(m + 1), x], x] /; FreeQ[{a, b, e, f, g}, x] && EqQ[a^2 - b^2, 0]
&& LtQ[m, -1] && GtQ[p, 1] && (GtQ[m, -2] || EqQ[2*m + p + 1, 0] || (EqQ[m, -2] && IntegerQ[p])) && NeQ[m + p,
 0] && IntegersQ[2*m, 2*p]

Rule 2759

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Simp[2*g*(g
*Cos[e + f*x])^(p - 1)*((a + b*Sin[e + f*x])^(m + 1)/(b*f*(2*m + p + 1))), x] + Dist[g^2*((p - 1)/(b^2*(2*m +
p + 1))), Int[(g*Cos[e + f*x])^(p - 2)*(a + b*Sin[e + f*x])^(m + 2), x], x] /; FreeQ[{a, b, e, f, g}, x] && Eq
Q[a^2 - b^2, 0] && LeQ[m, -2] && GtQ[p, 1] && NeQ[2*m + p + 1, 0] &&  !ILtQ[m + p + 1, 0] && IntegersQ[2*m, 2*
p]

Rule 2761

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)/((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Simp[g*((g*Cos[e
 + f*x])^(p - 1)/(b*f*(p - 1))), x] + Dist[g^2/a, Int[(g*Cos[e + f*x])^(p - 2), x], x] /; FreeQ[{a, b, e, f, g
}, x] && EqQ[a^2 - b^2, 0] && GtQ[p, 1] && IntegerQ[2*p]

Rule 2938

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.)
+ (f_.)*(x_)]), x_Symbol] :> Simp[(b*c - a*d)*(g*Cos[e + f*x])^(p + 1)*((a + b*Sin[e + f*x])^m/(a*f*g*(2*m + p
 + 1))), x] + Dist[(a*d*m + b*c*(m + p + 1))/(a*b*(2*m + p + 1)), Int[(g*Cos[e + f*x])^p*(a + b*Sin[e + f*x])^
(m + 1), x], x] /; FreeQ[{a, b, c, d, e, f, g, m, p}, x] && EqQ[a^2 - b^2, 0] && (LtQ[m, -1] || ILtQ[Simplify[
m + p], 0]) && NeQ[2*m + p + 1, 0]

Rule 3046

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Dist[a^m*c^m, Int[Cos[e + f*x]^(2*m)*(c + d*Sin[e + f*x])^(n - m)*(A + B
*Sin[e + f*x]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B, n}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && I
ntegerQ[m] &&  !(IntegerQ[n] && ((LtQ[m, 0] && GtQ[n, 0]) || LtQ[0, n, m] || LtQ[m, n, 0]))

Rubi steps \begin{align*} \text {integral}& = \left (a^3 c^3\right ) \int \frac {\cos ^6(e+f x) (A+B \sin (e+f x))}{(a+a \sin (e+f x))^5} \, dx \\ & = -\frac {a^3 (A-B) c^3 \cos ^7(e+f x)}{3 f (a+a \sin (e+f x))^5}-\frac {1}{3} \left (a^2 (2 A-5 B) c^3\right ) \int \frac {\cos ^6(e+f x)}{(a+a \sin (e+f x))^4} \, dx \\ & = -\frac {a^3 (A-B) c^3 \cos ^7(e+f x)}{3 f (a+a \sin (e+f x))^5}+\frac {2 a (2 A-5 B) c^3 \cos ^5(e+f x)}{3 f (a+a \sin (e+f x))^3}+\frac {1}{3} \left (5 (2 A-5 B) c^3\right ) \int \frac {\cos ^4(e+f x)}{(a+a \sin (e+f x))^2} \, dx \\ & = -\frac {a^3 (A-B) c^3 \cos ^7(e+f x)}{3 f (a+a \sin (e+f x))^5}+\frac {2 a (2 A-5 B) c^3 \cos ^5(e+f x)}{3 f (a+a \sin (e+f x))^3}+\frac {5 (2 A-5 B) c^3 \cos ^3(e+f x)}{6 f \left (a^2+a^2 \sin (e+f x)\right )}+\frac {\left (5 (2 A-5 B) c^3\right ) \int \frac {\cos ^2(e+f x)}{a+a \sin (e+f x)} \, dx}{2 a} \\ & = \frac {5 (2 A-5 B) c^3 \cos (e+f x)}{2 a^2 f}-\frac {a^3 (A-B) c^3 \cos ^7(e+f x)}{3 f (a+a \sin (e+f x))^5}+\frac {2 a (2 A-5 B) c^3 \cos ^5(e+f x)}{3 f (a+a \sin (e+f x))^3}+\frac {5 (2 A-5 B) c^3 \cos ^3(e+f x)}{6 f \left (a^2+a^2 \sin (e+f x)\right )}+\frac {\left (5 (2 A-5 B) c^3\right ) \int 1 \, dx}{2 a^2} \\ & = \frac {5 (2 A-5 B) c^3 x}{2 a^2}+\frac {5 (2 A-5 B) c^3 \cos (e+f x)}{2 a^2 f}-\frac {a^3 (A-B) c^3 \cos ^7(e+f x)}{3 f (a+a \sin (e+f x))^5}+\frac {2 a (2 A-5 B) c^3 \cos ^5(e+f x)}{3 f (a+a \sin (e+f x))^3}+\frac {5 (2 A-5 B) c^3 \cos ^3(e+f x)}{6 f \left (a^2+a^2 \sin (e+f x)\right )} \\ \end{align*}

Mathematica [A] (verified)

Time = 11.30 (sec) , antiderivative size = 274, normalized size of antiderivative = 1.69 \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=\frac {\left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right ) (c-c \sin (e+f x))^3 \left (64 (A-B) \sin \left (\frac {1}{2} (e+f x)\right )-32 (A-B) \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )-32 (7 A-13 B) \sin \left (\frac {1}{2} (e+f x)\right ) \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^2+30 (2 A-5 B) (e+f x) \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^3+12 (A-5 B) \cos (e+f x) \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^3+3 B \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^3 \sin (2 (e+f x))\right )}{12 a^2 f \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^6 (1+\sin (e+f x))^2} \]

[In]

Integrate[((A + B*Sin[e + f*x])*(c - c*Sin[e + f*x])^3)/(a + a*Sin[e + f*x])^2,x]

[Out]

((Cos[(e + f*x)/2] + Sin[(e + f*x)/2])*(c - c*Sin[e + f*x])^3*(64*(A - B)*Sin[(e + f*x)/2] - 32*(A - B)*(Cos[(
e + f*x)/2] + Sin[(e + f*x)/2]) - 32*(7*A - 13*B)*Sin[(e + f*x)/2]*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^2 + 3
0*(2*A - 5*B)*(e + f*x)*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^3 + 12*(A - 5*B)*Cos[e + f*x]*(Cos[(e + f*x)/2]
+ Sin[(e + f*x)/2])^3 + 3*B*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^3*Sin[2*(e + f*x)]))/(12*a^2*f*(Cos[(e + f*x
)/2] - Sin[(e + f*x)/2])^6*(1 + Sin[e + f*x])^2)

Maple [A] (verified)

Time = 0.94 (sec) , antiderivative size = 162, normalized size of antiderivative = 1.00

method result size
derivativedivides \(\frac {2 c^{3} \left (\frac {-\frac {B \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2}+\left (A -5 B \right ) \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\frac {B \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{2}+A -5 B}{\left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}+\frac {5 \left (2 A -5 B \right ) \arctan \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2}-\frac {-16 A +16 B}{2 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{2}}-\frac {-4 A +12 B}{\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1}-\frac {16 A -16 B}{3 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}}\right )}{f \,a^{2}}\) \(162\)
default \(\frac {2 c^{3} \left (\frac {-\frac {B \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2}+\left (A -5 B \right ) \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\frac {B \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{2}+A -5 B}{\left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}+\frac {5 \left (2 A -5 B \right ) \arctan \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2}-\frac {-16 A +16 B}{2 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{2}}-\frac {-4 A +12 B}{\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1}-\frac {16 A -16 B}{3 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}}\right )}{f \,a^{2}}\) \(162\)
parallelrisch \(-\frac {c^{3} \left (\left (-30 f x A +75 f x B -56 A +\frac {607}{4} B \right ) \cos \left (\frac {f x}{2}+\frac {e}{2}\right )+\left (10 f x A -25 f x B -\frac {19}{3} A +\frac {115}{12} B \right ) \cos \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )+\left (-30 f x A +75 f x B -36 A +\frac {337}{4} B \right ) \sin \left (\frac {f x}{2}+\frac {e}{2}\right )+\left (-10 f x A +25 f x B -37 A +\frac {353}{4} B \right ) \sin \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )+\left (A -\frac {17 B}{4}\right ) \cos \left (\frac {5 f x}{2}+\frac {5 e}{2}\right )+\left (-A +\frac {17 B}{4}\right ) \sin \left (\frac {5 f x}{2}+\frac {5 e}{2}\right )+\frac {B \left (\cos \left (\frac {7 f x}{2}+\frac {7 e}{2}\right )+\sin \left (\frac {7 f x}{2}+\frac {7 e}{2}\right )\right )}{4}\right )}{2 f \,a^{2} \left (3 \cos \left (\frac {f x}{2}+\frac {e}{2}\right )-\cos \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )+\sin \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )+3 \sin \left (\frac {f x}{2}+\frac {e}{2}\right )\right )}\) \(220\)
risch \(\frac {5 c^{3} x A}{a^{2}}-\frac {25 c^{3} x B}{2 a^{2}}-\frac {i B \,c^{3} {\mathrm e}^{2 i \left (f x +e \right )}}{8 a^{2} f}+\frac {c^{3} {\mathrm e}^{i \left (f x +e \right )} A}{2 a^{2} f}-\frac {5 c^{3} {\mathrm e}^{i \left (f x +e \right )} B}{2 a^{2} f}+\frac {c^{3} {\mathrm e}^{-i \left (f x +e \right )} A}{2 a^{2} f}-\frac {5 c^{3} {\mathrm e}^{-i \left (f x +e \right )} B}{2 a^{2} f}+\frac {i B \,c^{3} {\mathrm e}^{-2 i \left (f x +e \right )}}{8 a^{2} f}+\frac {32 i A \,c^{3} {\mathrm e}^{i \left (f x +e \right )}+24 A \,c^{3} {\mathrm e}^{2 i \left (f x +e \right )}-64 i B \,c^{3} {\mathrm e}^{i \left (f x +e \right )}-40 B \,c^{3} {\mathrm e}^{2 i \left (f x +e \right )}-\frac {56 A \,c^{3}}{3}+\frac {104 B \,c^{3}}{3}}{f \,a^{2} \left ({\mathrm e}^{i \left (f x +e \right )}+i\right )^{3}}\) \(247\)
norman \(\frac {\frac {\left (8 A \,c^{3}-25 B \,c^{3}\right ) \left (\tan ^{10}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a f}+\frac {\left (34 A \,c^{3}-77 B \,c^{3}\right ) \left (\tan ^{9}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a f}+\frac {\left (38 A \,c^{3}-93 B \,c^{3}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{a f}+\frac {\left (148 A \,c^{3}-352 B \,c^{3}\right ) \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a f}+\frac {46 A \,c^{3}-118 B \,c^{3}}{3 a f}+\frac {5 c^{3} \left (2 A -5 B \right ) x}{2 a}+\frac {2 \left (68 A \,c^{3}-194 B \,c^{3}\right ) \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a f}+\frac {2 \left (70 A \,c^{3}-160 B \,c^{3}\right ) \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a f}+\frac {2 \left (108 A \,c^{3}-251 B \,c^{3}\right ) \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a f}+\frac {\left (154 A \,c^{3}-478 B \,c^{3}\right ) \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 a f}+\frac {2 \left (182 A \,c^{3}-545 B \,c^{3}\right ) \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 a f}+\frac {\left (220 A \,c^{3}-595 B \,c^{3}\right ) \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 a f}+\frac {15 c^{3} \left (2 A -5 B \right ) x \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{2 a}+\frac {35 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2 a}+\frac {65 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2 a}+\frac {45 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a}+\frac {55 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a}+\frac {55 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a}+\frac {45 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a}+\frac {65 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2 a}+\frac {35 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{9}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2 a}+\frac {15 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{10}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2 a}+\frac {5 c^{3} \left (2 A -5 B \right ) x \left (\tan ^{11}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2 a}}{\left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{4} a \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}}\) \(680\)

[In]

int((A+B*sin(f*x+e))*(c-c*sin(f*x+e))^3/(a+a*sin(f*x+e))^2,x,method=_RETURNVERBOSE)

[Out]

2/f*c^3/a^2*((-1/2*B*tan(1/2*f*x+1/2*e)^3+(A-5*B)*tan(1/2*f*x+1/2*e)^2+1/2*B*tan(1/2*f*x+1/2*e)+A-5*B)/(1+tan(
1/2*f*x+1/2*e)^2)^2+5/2*(2*A-5*B)*arctan(tan(1/2*f*x+1/2*e))-1/2*(-16*A+16*B)/(tan(1/2*f*x+1/2*e)+1)^2-(-4*A+1
2*B)/(tan(1/2*f*x+1/2*e)+1)-1/3*(16*A-16*B)/(tan(1/2*f*x+1/2*e)+1)^3)

Fricas [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 291, normalized size of antiderivative = 1.80 \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=\frac {3 \, B c^{3} \cos \left (f x + e\right )^{4} + 6 \, {\left (A - 4 \, B\right )} c^{3} \cos \left (f x + e\right )^{3} - 30 \, {\left (2 \, A - 5 \, B\right )} c^{3} f x + 16 \, {\left (A - B\right )} c^{3} + {\left (15 \, {\left (2 \, A - 5 \, B\right )} c^{3} f x - {\left (62 \, A - 131 \, B\right )} c^{3}\right )} \cos \left (f x + e\right )^{2} - {\left (15 \, {\left (2 \, A - 5 \, B\right )} c^{3} f x + 2 \, {\left (26 \, A - 71 \, B\right )} c^{3}\right )} \cos \left (f x + e\right ) + {\left (3 \, B c^{3} \cos \left (f x + e\right )^{3} - 30 \, {\left (2 \, A - 5 \, B\right )} c^{3} f x - 3 \, {\left (2 \, A - 9 \, B\right )} c^{3} \cos \left (f x + e\right )^{2} - 16 \, {\left (A - B\right )} c^{3} - {\left (15 \, {\left (2 \, A - 5 \, B\right )} c^{3} f x + 2 \, {\left (34 \, A - 79 \, B\right )} c^{3}\right )} \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )}{6 \, {\left (a^{2} f \cos \left (f x + e\right )^{2} - a^{2} f \cos \left (f x + e\right ) - 2 \, a^{2} f - {\left (a^{2} f \cos \left (f x + e\right ) + 2 \, a^{2} f\right )} \sin \left (f x + e\right )\right )}} \]

[In]

integrate((A+B*sin(f*x+e))*(c-c*sin(f*x+e))^3/(a+a*sin(f*x+e))^2,x, algorithm="fricas")

[Out]

1/6*(3*B*c^3*cos(f*x + e)^4 + 6*(A - 4*B)*c^3*cos(f*x + e)^3 - 30*(2*A - 5*B)*c^3*f*x + 16*(A - B)*c^3 + (15*(
2*A - 5*B)*c^3*f*x - (62*A - 131*B)*c^3)*cos(f*x + e)^2 - (15*(2*A - 5*B)*c^3*f*x + 2*(26*A - 71*B)*c^3)*cos(f
*x + e) + (3*B*c^3*cos(f*x + e)^3 - 30*(2*A - 5*B)*c^3*f*x - 3*(2*A - 9*B)*c^3*cos(f*x + e)^2 - 16*(A - B)*c^3
 - (15*(2*A - 5*B)*c^3*f*x + 2*(34*A - 79*B)*c^3)*cos(f*x + e))*sin(f*x + e))/(a^2*f*cos(f*x + e)^2 - a^2*f*co
s(f*x + e) - 2*a^2*f - (a^2*f*cos(f*x + e) + 2*a^2*f)*sin(f*x + e))

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4665 vs. \(2 (148) = 296\).

Time = 7.42 (sec) , antiderivative size = 4665, normalized size of antiderivative = 28.80 \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=\text {Too large to display} \]

[In]

integrate((A+B*sin(f*x+e))*(c-c*sin(f*x+e))**3/(a+a*sin(f*x+e))**2,x)

[Out]

Piecewise((30*A*c**3*f*x*tan(e/2 + f*x/2)**7/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 3
0*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e
/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 90*A*c**3*f*x*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2
+ f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 +
42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 150*A
*c**3*f*x*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/
2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2
+ 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 210*A*c**3*f*x*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 +
18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(
e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 210*A*c**3*f*x*tan(
e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5
+ 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*ta
n(e/2 + f*x/2) + 6*a**2*f) + 150*A*c**3*f*x*tan(e/2 + f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(
e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**
3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 90*A*c**3*f*x*tan(e/2 + f*x/2)/(6
*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/
2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6
*a**2*f) + 30*A*c**3*f*x/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f
*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*
a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 48*A*c**3*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*
tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/
2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 204*A*c**3*tan(e/2 + f*x/2)**
5/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*ta
n(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2)
 + 6*a**2*f) + 212*A*c**3*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 +
30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(
e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 432*A*c**3*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 +
f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42
*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 256*A*c
**3*tan(e/2 + f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*
x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a
**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 228*A*c**3*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan
(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)*
*3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 92*A*c**3/(6*a**2*f*tan(e/2 + f*
x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a
**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 75*B*c**3
*f*x*tan(e/2 + f*x/2)**7/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f
*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*
a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 225*B*c**3*f*x*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a*
*2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 +
 f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 375*B*c**3*f*x*tan(e/2 +
 f*x/2)**5/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*
a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2
 + f*x/2) + 6*a**2*f) - 525*B*c**3*f*x*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 +
 f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 3
0*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 525*B*c**3*f*x*tan(e/2 + f*x/2)**3/(6*
a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2
 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*
a**2*f) - 375*B*c**3*f*x*tan(e/2 + f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 3
0*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e
/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 225*B*c**3*f*x*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 +
f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42
*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 75*B*c*
*3*f*x/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2
*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f
*x/2) + 6*a**2*f) - 150*B*c**3*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)*
*6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f
*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 462*B*c**3*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e
/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4
 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 65
6*B*c**3*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2
 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 +
 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 996*B*c**3*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a*
*2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 +
 f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 718*B*c**3*tan(e/2 + f*x
/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2
*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f
*x/2) + 6*a**2*f) - 558*B*c**3*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6
+ 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*ta
n(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 236*B*c**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**
2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 +
f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f), Ne(f, 0)), (x*(A + B*sin(e
))*(-c*sin(e) + c)**3/(a*sin(e) + a)**2, True))

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1378 vs. \(2 (152) = 304\).

Time = 0.33 (sec) , antiderivative size = 1378, normalized size of antiderivative = 8.51 \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=\text {Too large to display} \]

[In]

integrate((A+B*sin(f*x+e))*(c-c*sin(f*x+e))^3/(a+a*sin(f*x+e))^2,x, algorithm="maxima")

[Out]

-1/3*(B*c^3*((75*sin(f*x + e)/(cos(f*x + e) + 1) + 97*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 126*sin(f*x + e)^3
/(cos(f*x + e) + 1)^3 + 98*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 63*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 21*s
in(f*x + e)^6/(cos(f*x + e) + 1)^6 + 32)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 5*a^2*sin(f*x + e)^2/(
cos(f*x + e) + 1)^2 + 7*a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 7*a^2*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 +
5*a^2*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 3*a^2*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 + a^2*sin(f*x + e)^7/(co
s(f*x + e) + 1)^7) + 21*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/a^2) - 4*A*c^3*((12*sin(f*x + e)/(cos(f*x + e)
 + 1) + 11*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 9*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 3*sin(f*x + e)^4/(cos
(f*x + e) + 1)^4 + 5)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 4*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1)^2
 + 4*a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 3*a^2*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + a^2*sin(f*x + e)^5/
(cos(f*x + e) + 1)^5) + 3*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/a^2) + 12*B*c^3*((12*sin(f*x + e)/(cos(f*x +
 e) + 1) + 11*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 9*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 3*sin(f*x + e)^4/(
cos(f*x + e) + 1)^4 + 5)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 4*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1
)^2 + 4*a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 3*a^2*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + a^2*sin(f*x + e)
^5/(cos(f*x + e) + 1)^5) + 3*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/a^2) - 6*A*c^3*((9*sin(f*x + e)/(cos(f*x
+ e) + 1) + 3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 4)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 3*a^2*si
n(f*x + e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3) + 3*arctan(sin(f*x + e)/(cos(f*x
+ e) + 1))/a^2) + 6*B*c^3*((9*sin(f*x + e)/(cos(f*x + e) + 1) + 3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 4)/(a^
2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 3*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f*x + e)^3/(co
s(f*x + e) + 1)^3) + 3*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/a^2) + 2*A*c^3*(3*sin(f*x + e)/(cos(f*x + e) +
1) + 3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 2)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 3*a^2*sin(f*x +
 e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3) - 6*A*c^3*(3*sin(f*x + e)/(cos(f*x + e)
+ 1) + 1)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 3*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f
*x + e)^3/(cos(f*x + e) + 1)^3) + 2*B*c^3*(3*sin(f*x + e)/(cos(f*x + e) + 1) + 1)/(a^2 + 3*a^2*sin(f*x + e)/(c
os(f*x + e) + 1) + 3*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3))/f

Giac [A] (verification not implemented)

none

Time = 0.31 (sec) , antiderivative size = 222, normalized size of antiderivative = 1.37 \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=\frac {\frac {15 \, {\left (2 \, A c^{3} - 5 \, B c^{3}\right )} {\left (f x + e\right )}}{a^{2}} - \frac {6 \, {\left (B c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} - 2 \, A c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 10 \, B c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - B c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 2 \, A c^{3} + 10 \, B c^{3}\right )}}{{\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 1\right )}^{2} a^{2}} + \frac {16 \, {\left (3 \, A c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - 9 \, B c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 12 \, A c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 24 \, B c^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 5 \, A c^{3} - 11 \, B c^{3}\right )}}{a^{2} {\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 1\right )}^{3}}}{6 \, f} \]

[In]

integrate((A+B*sin(f*x+e))*(c-c*sin(f*x+e))^3/(a+a*sin(f*x+e))^2,x, algorithm="giac")

[Out]

1/6*(15*(2*A*c^3 - 5*B*c^3)*(f*x + e)/a^2 - 6*(B*c^3*tan(1/2*f*x + 1/2*e)^3 - 2*A*c^3*tan(1/2*f*x + 1/2*e)^2 +
 10*B*c^3*tan(1/2*f*x + 1/2*e)^2 - B*c^3*tan(1/2*f*x + 1/2*e) - 2*A*c^3 + 10*B*c^3)/((tan(1/2*f*x + 1/2*e)^2 +
 1)^2*a^2) + 16*(3*A*c^3*tan(1/2*f*x + 1/2*e)^2 - 9*B*c^3*tan(1/2*f*x + 1/2*e)^2 + 12*A*c^3*tan(1/2*f*x + 1/2*
e) - 24*B*c^3*tan(1/2*f*x + 1/2*e) + 5*A*c^3 - 11*B*c^3)/(a^2*(tan(1/2*f*x + 1/2*e) + 1)^3))/f

Mupad [B] (verification not implemented)

Time = 14.79 (sec) , antiderivative size = 336, normalized size of antiderivative = 2.07 \[ \int \frac {(A+B \sin (e+f x)) (c-c \sin (e+f x))^3}{(a+a \sin (e+f x))^2} \, dx=\frac {\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (38\,A\,c^3-93\,B\,c^3\right )+\frac {46\,A\,c^3}{3}-\frac {118\,B\,c^3}{3}+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^6\,\left (8\,A\,c^3-25\,B\,c^3\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^5\,\left (34\,A\,c^3-77\,B\,c^3\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^3\,\left (72\,A\,c^3-166\,B\,c^3\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^4\,\left (\frac {106\,A\,c^3}{3}-\frac {328\,B\,c^3}{3}\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2\,\left (\frac {128\,A\,c^3}{3}-\frac {359\,B\,c^3}{3}\right )}{f\,\left (a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^7+3\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^6+5\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^5+7\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^4+7\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^3+5\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2+3\,a^2\,\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )+a^2\right )}+\frac {5\,c^3\,\mathrm {atan}\left (\frac {5\,c^3\,\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (2\,A-5\,B\right )}{10\,A\,c^3-25\,B\,c^3}\right )\,\left (2\,A-5\,B\right )}{a^2\,f} \]

[In]

int(((A + B*sin(e + f*x))*(c - c*sin(e + f*x))^3)/(a + a*sin(e + f*x))^2,x)

[Out]

(tan(e/2 + (f*x)/2)*(38*A*c^3 - 93*B*c^3) + (46*A*c^3)/3 - (118*B*c^3)/3 + tan(e/2 + (f*x)/2)^6*(8*A*c^3 - 25*
B*c^3) + tan(e/2 + (f*x)/2)^5*(34*A*c^3 - 77*B*c^3) + tan(e/2 + (f*x)/2)^3*(72*A*c^3 - 166*B*c^3) + tan(e/2 +
(f*x)/2)^4*((106*A*c^3)/3 - (328*B*c^3)/3) + tan(e/2 + (f*x)/2)^2*((128*A*c^3)/3 - (359*B*c^3)/3))/(f*(5*a^2*t
an(e/2 + (f*x)/2)^2 + 7*a^2*tan(e/2 + (f*x)/2)^3 + 7*a^2*tan(e/2 + (f*x)/2)^4 + 5*a^2*tan(e/2 + (f*x)/2)^5 + 3
*a^2*tan(e/2 + (f*x)/2)^6 + a^2*tan(e/2 + (f*x)/2)^7 + a^2 + 3*a^2*tan(e/2 + (f*x)/2))) + (5*c^3*atan((5*c^3*t
an(e/2 + (f*x)/2)*(2*A - 5*B))/(10*A*c^3 - 25*B*c^3))*(2*A - 5*B))/(a^2*f)